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Abstract 

We have searched a sample of 9.6 million BB events for the lepton-fiavor-violating leptonic 
B decays, B° im^t^ and B^ — > e='=r^. The r-lepton was detected through the decay modes 

T ivv, where t = e,[i. There is no indication of a signal, and we obtain the 90% confidence level 
upper limits B{B^ h^t^) < 3.8 x 10"^ and ^ e^r^) < 1.1 x 10"^. 



*0n leave of absence from University of Chicago. 
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We report results of a search for two lepton-flavor-violating leptonic decays of B mesons: 
/i^r^ and B^ — > e'^r^. These modes are forbidden in the conventional Standard 
Model by the lepton-flavor conservation law. However, they are predicted to occur in many 
theories "beyond the standard model," for example multi-Higgs-boson extensions, theories 
with leptoquarks, supersymmetric models without R-parity, and Higgs-mediated decay in su- 
persymmetric seesaw models fl'|. The recent discovery of neutrino oscillation, while not lead- 
ing to predictions of observable rates for lepton-flavor-violating decays, nonetheless heightens 
interest in them The decays we searched for involve both third generation quarks and 
third generation leptons. Decays of this variety have been less extensively searched for than 
those involving only first or second generation quarks or leptons. Discovery of such decays 
at levels of our sensitivity would be clear evidence of physics beyond the Standard Model. 
Currently the best limits on the branching fractions are B{B^ — > /x^r^) < 8.3 x 10~^, and 
B{B° e^T^) < 5.3 X 10-^ at 90% confidence level 0. 

The data used in this analysis were taken with the CLEO detector 4] at the Cornell 
Electron Storage Ring (CESR), a symmetric e~^e~ collider operating in the T(4S) resonance 
region. The data sample consists of 9.2 fb~^ at the resonance, corresponding to 9.6 million 
BB events, and 4.5 fb~^ at a center-of-mass energy 60 MeV below the resonance. The sample 
below the resonance provides information on the background from continuum processes 
e~^e~ ^ qq, q = u,d,s,c, and from two-photon fusion processes e'^e" e"'"e~7*7*, 7*7* — >■ 
X (7* a virtual photon). We scale the off-resonance yields by 1.99, the luminosity ratio 
divided by the cm. energy-squared ratio, and subtract them from on-resonance yields to 
obtain BB yields. 

Summing over fi{e)~^T~ and ^{e)~T~^, we search for 5° fi{e)^T^ with the r lepton 
detected via the r — euu and r fiuu decay modes. In this article, i denotes the primary 
lepton from the signal B and i' denotes the secondary lepton from r. {i, i') denotes B^ — *■ 
£^r^, T £Vz/. We have four modes to analyze; {fi, e), {fi, /x), (e, e), and (e, fi). 

Muons are identified by their ability to penetrate the iron return yoke of the magnet: at 
least five (three) interaction lengths of material for the primary (secondary) muon. Electrons 
are identified by shower energy to momentum ratio [E/P), track-cluster matching, dE/dx, 
and shower shape. 

In the rest frame of the signal B, the primary lepton is monoenergetic, with momentum 
2.34 GeV/c. In the lab frame (the T(4S) rest frame), this is smeared, and ranges from 2.2 
to 2.5 GeV/c. We require that the primary lepton candidate have momentum in that range. 
We require that the secondary lepton, from r, be greater than 0.6 (1.0) GeV/c for e (/x). 
We "measure" the 4-momentum of the neutrino pair as the missing visible 4-momentum in 
the event : Ej^p = 2Ebeam — ^Ei, P^p = —T,Pi, where sums are over all observed (charged 
and neutral) particles. 

We define two neural net variables. NN^j^ is a neural net variable used to sup- 
press backgrounds from BB decays. We calculate three inputs: beam-constrained mass 

\jEtam - Pcand^ = E^and " ^feeam, whcrc P^and (Ecand) IS the momeutum (energy) of the 
B candidate, and cos 613 (the cosine of the angle between the momenta of primary lepton 
and B candidate). We feed them into a neural net and train it with signal and BB Monte 
Carlo simulations for each mode. NN^ont is a neural net variable to suppress backgrounds 
from continuum. We calculate five inputs: R2 (the ratio of second and zeroth Fox- Wolfram 
moments jsj of the event), S (the sphericity), thrust of the event, cos 9u (the cosine of the 
angle between the pi — pi' and the thrust axis of the rest of the event), and cos 9^^-^^^+^^, 
(the cosine of the angle between the momenta of neutrino pair and lepton pair), then feed 
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them into a neural net and train it with signal and continuum Monte Carlo for each mode. 
The nominal neural net range is from 0.0 to 1.0. We cut in the 2D space defined by NN^q 
and NN^onu requiring A^A^^^ > NN^j^^, NN^ont > NN^^^^ and also 

We define two r-mass variables. The first is the conventionally defined invariant mass 
of the reconstructed r, M^vp = \/ {Eg' + E^pY — [pii + p^p)^. The second r-mass variable 
makes use of the fact that, with perfect measurements of all quantities, AE = 0, and hence 
we can use Ebeam - Eg - Ei> for E^^p, yielding Mi>^p^AE=o = \/ {Ebeam - EeY - {pe + Puu^. 
We further define AM^ = Mp^p — Mr and AM^^a£;=o = ^r;/!5,AE=o — M^, where Mr is the 
nominal r mass, 1777 MeV. 

By examining the angular distribution of electrons, positrons, and missing momentum, 
in off-resonance data in a \pe\ sideband region (2.0 < \pe\ < 2.2 GeV/c and 2.5 < \pe\ < 2.7 
GeV/c), we see clear evidence of the two-photon-fusion process. There are sharp peaks in 
the forward directions for electrons and positrons ( "forward" being the direction of the beam 
particle of the same charge). Also, the missing momentum peaks sharply in the opposite 
direction from a detected e"*" or e~, indicating an e~ or e"*" lost down the beam pipe. By 
eliminating events with \ cos9miss\ > 0.90 (0.95 for (/i, /i)), we considerably reduce this 
background. 

We compare Monte Carlo samples with data using the \pe\ sideband region defined above. 
In Fig. HJ we show distributions in NNcont, NN^^, AMr, and AMt-^ae=o, for off-resonance- 
subtracted on-resonance data and absolutely normalized BB Monte Carlo. Agreement is 
good. In Fig. |21 we show distributions for the same variables for off-resonance data and 
absolutely normalized continuum {e~^e~ ^ qq, q = u, d, s, c) Monte Carlo. We have not 
included a Monte Carlo for the inclusive multi-hadronic two-photon-fusion process, lacking 
a trustworthy simulation of this process. For (/x, e), shown in Fig. |2l agreement is good, 
indicating that the remaining contribution from two-photon fusion is small. For (e, e), not 
shown, data exceeds continuum Monte Carlo, indicating a sizeable remaining contribution 
from two-photon fusion. 

We measure the ratio of data to Monte Carlo yields in the \pe\ sideband region, denoting 
hy TZ^B the ratio of off- resonance-subtracted on-resonance data to BB Monte Carlo, and by 
TZcont the ratio of off-resonance data to continuum Monte Carlo. 7^'s are measured with loose 
selection criteria applied: NNcont > 0.5, NN^j^ > 0.5, \ cos9miss\ < 0.9 (0.95 for (/i, fi)), 
and |AM^| < 2.0 GeV. Values so obtained are given in Table E One sees that TZbb differs 
little from 1.0, while TZcont is less well behaved, particularly for the (e, e) mode, indicating 
that the two-photon fusion background is present. To estimate backgrounds in the signal 
\pi\ region correctly, we scale the signal- region yields from Monte Carlo by TZ. Because we 
perform a direct subtraction of off-resonance data, the accuracy of the continuum background 
prediction is not critical for our results. 

We optimized our selection criteria on NN^s, NNcont, AMr, and AM^^a_b=o to obtain 
the best upper limit when the true branching fraction is zero. This optimization procedure 
made use of signal and background Monte Carlo samples, and scaled the background samples 
by TZ^B or TZcont as described above. The optimized selection criteria, found separately for 
each mode, are shown in Table lU 

The number of events satisfying all selection criteria is shown, for each mode, in Table HI 
along with the background estimate. We find 15 {fi, e) candidates, with 23.7 expected from 
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FIG. 1: (a) NNconti (b) NN^j^, (c) AMr, and (d) AMr^A.E=o distributions, for the comparison 
of BB Monte Carlo (histogram) vs off-resonance-subtracted on-resonance data (points) in the p£ 
sideband region, for the (/i, e) mode. 




FIG. 2: (a) NNcont, (b) NN^g, (c) AM,-, and (d) AMt-^a_b=o distributions, for the comparison 
of continuum Monte Carlo (e^e~ qq, q = u,d, s, c: histogram) vs off-resonance data (points) in 
the Pi sideband region, for the (/i, e) mode. 



background; we find 4 (/i, /i) candidates, with 9.0 expected from background; we find 14 
(e, e) candidates, with 11.6 expected from background; we find 6 (e, /i) candidates, with 
5.1 expected from background. Thus there are a total of 39 events with 49.4 expected from 
background. The probability that a true mean of 49.4 will give rise to a yield of 39 or more 
events is 93%. With no indication of signal, we obtain the branching fraction upper limits. 

We calculate upper limits at 90% confidence level. There is some probability of observing 
the off-resonance-subtracted on-resonance yield that we do observe, or less, if the branch- 
ing fractions for /x^r^ and — > e^r^ are zero. We take the 90% confidence level 
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TABLE I: The rows of NN^^, NN^^^^, AM^ (GeV), and AMr,AE=o (GeV) indicate the optimized 
selection criteria for each mode. TZbb T^cont are measured ratios between data yields and 
Monte Carlo with loose selection criteria in the \pe\ sideband region. Non{Noff) is the number of 
observed events satisfying the optimized selection criteria in the signal region of from on(off)- 
resonance data samples. Nobs is the number of observed events from off-resonance-subtracted 
on-resonance data, Nqn—^-^^Noff- is the BB background estimate from Monte Carlo, 

scaled by TZ^g] in the absence of signal it should be comparable with N^bs- O-^Ncont is (1/1.99) 
times the continuum background estimate from Monte Carlo, scaled by TZcont, which should be 
comparable with Nqff- e is the signal detection efficiency including r decay branching fraction. 
BR UL is the branching ratio upper limit at 90% confidence level with systematic error considered. 





(/^, e) 


IJ) 


(e, e) 


(e,)u) 


-2.0<AM^ 
-2.0<AM^,Ai<;=o 


0.725 
0.700 
<2.00 
<0.25 


0.875 
0.775 
<1.40 
<0.25 


0.675 
0.700 
<1.50 
<0.30 


0.825 
0.475 
<1.40 
<0.25 




1.21±0.06 
1.03±0.27 


1.06±0.07 
1.52±0.32 


1.04±0.07 
4.30±0.95 


0.94±0.07 
0.64±0.38 


NON 

Nqff 

Nobs 


19 
2 

15.0 ± 5.2 


10 
3 

4.0±4.7 


28 
7 

14.0 ±7.5 


6 


6.0 ± 2.4 


Nbb 
0.57V^„t 


23.7±2.7 
1.8±0.6 


9.0±1.4 
0.4±0.2 


11.6 ±1.4 
3.1±1.0 


5.1 ±0.8 
0.5±0.3 


6(%) 


1.57 


0.63 


0.96 


0.58 


BR UL{IQ-^) 


0.55 


0.87 


1.64 


1.46 



upper limit to be that value of the branching fraction which reduces the above-mentioned 

probability by a factor of 10. The ingredients needed for the calculation are: (1) the ob- 
served ofT-resonance-subtracted on-resonance yield; (2) the true mean for the background 
contribution from BB processes, and (3) less critically, the true mean of the background 
contribution from non-resonance processes. To allow for the uncertainty in the background 
estimates, we changed TZ^b ^-nd TZcont in the unfavorable directions by Icr, i.e. —la for TZbb 
and +1(7 for IZcont- 

We use Monte Carlo simulation to determine the efficiency for detecting the signal modes. 
The decays B — > ^{e)^T^ are generated with the r lepton unpolarized. For a r lepton 
polarization as given by — ^, the secondary lepton is boosted (has its average lab energy 
increased), which in turn increases the efficiency. For the opposite polarization, as given 
\yy V + A.I the secondary lepton is deboosted, and the efficiency is lowered. The fractional 
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changes in efficiency, averaged over the four modes, are +11% for V — A, —8% for V + A. 
Our upper hmits are quoted for unpolarized r's. 

Systematic errors are of two varieties - those on the estimate of signal detection efficien- 
cies, and those on the estimate of backgrounds. The dominant contributors to the former 
are lepton identification efficiency uncertainties (contributing ±3.5% per lepton, relative, in 
the efficiency), and missing-four- vector-simulation uncertainties (±5.4%), giving a relative 
uncertainty in the overall efficiency of ±7.4% for (e, /i) and (/i, e), and ±8.9% for (e, e) 
and (/i, fi). The background uncertainties are handled by varying the TZbb Tlcont as 
mentioned above. The errors shown on the backgrounds in Table H] include statistical and 
systematic errors. 

There is no universally agreed-upon procedure for including systematic errors in upper- 
limit estimates. We conservatively vary the background by 1.0 standard deviations, and 
decrease the efficiency by 1.0 standard deviations for each mode and the results are as 
shown in Table EJ 

To combine the results from two leptonic modes, r — »• ew and r — > /ii/z/, we simply add 
the yields, add the backgrounds, and add the efficiencies. 
In this way we obtain our final results: 

B{B ^^T^) < 3.8 X 10"^ , 

B{B e^T^) < 1.1 X 10"^ , 

both at 90% confidence level. These results are significant improvements over previously 
pubhshed limits j^. 

In summary, we have searched for the decays B — >■ /i^r"^ and B e^r^. We find no 
indication of a signal, and obtain upper limits on the branching fractions. 

We gratefully acknowledge the effort of the CESR staff in providing us with excellent 
luminosity and running conditions. This work was supported by the National Science Foun- 
dation, the U.S. Department of Energy, the Research Corporation, and the Texas Advanced 
Research Program. 
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